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Spiral waves anchored to obstacles in cardiac tissues may cause lethal arrhythmia. To unpin these anchored 
spirals, comparing to high-voltage side-effect traditional therapies, wave emission from heterogeneities 
(WEH) induced by the uniform electric field (UEF) has provided a low-voltage alternative. Here we provide 
a new approach using WEH induced by the circularly polarized electric field (CPEF), which has higher 
success rate and larger application scope than UEF, even with a lower voltage. And we also study the 
distribution of the membrane potential near an obstacle induced by CPEF to analyze its mechanism of 
unpinning. We hope this promising approach may provide a better alternative to terminate arrhythmia. 



Spirals, also known as rotors^ or vortices^, occur in various excitable systems, including chemical media^"^, 
aggregations of Dictyostelium discoideum amoebae^, and cardiac tissues^. In hearts, spirals and subsequent 
turbulences may cause lethal arrhythmia^"^^. Better than traditional therapies^^ a new approach using wave 
emission from heterogeneities (WEH) or far-field stimulation provides a promising alternative to terminate 
arrhythmia^'*"^^. This approach is based on the fact that, by the application of an external electric field to a whole 
piece of tissue, de-polarizations and hyper-polarizations (so-called Weidmann zones^^) could be induced near 
obstacles (conductivity heterogeneities). These obstacles correspond to blood vessels, ischemic regions, and 
smaller- scale discontinuities. If the electric strength exceeds some threshold, obstacles can act as virtual electrodes 
or second sources^^"^^ 

The life-saving motivation to terminate arrhythmia has sparked many discussions about the mechanism of 
WEH^^"^^. However, previous works focus on WEH in response to the uniform electric field (UEF), which is 
realized by applying DC pulses onto field electrodes^^"^^. Recently, the circularly polarized electric field (CPEF) 
has shown its unique ability to control spirals and turbulences^^'^^, and has been verified in the Belousov- 
Zhabotinsky reaction by applying two ACs onto two pairs of field electrodes perpendicular to each other^^. 

In this paper, we study the mechanism of WEH induced by CPEF, and find that its ability to unpin anchored 
spirals, which is an important step in terminating arrhythmia, has advantages over UEF, such as lower voltage, 
higher success rate and larger application scope. Therefore, as a lower-voltage higher- efficiency approach, CPEF 
is more applicable in terminating arrhythmia. 

In the following, without loss of generality, we use a counter-clockwise rotating CPEF to unpin anchored 
spirals, which can be expressed as E = (£^, Ey), where = Eq cos{cOet + (/>e), Ey = Eq cos{cOet + (/>e + 37r/2) and Eq, 
cOe, are its strength, angular frequency and initial phase relative to x axis. In mono-domain models, the general 
effect of an external electric field on an obstacle can be expressed as an additional no -flux boundary condition^^'^^: 
n • V(y + E-r) = 0, where n is the normal vector to the obstacle boundary. Vis the membrane potential, E is the 
external electric field, and r is a point on the boundary. Therefore, in the presence of a circular obstacle (radius R) 
influenced by CPEF, the boundary condition can be described in the polar coordinate (p, 0) as 



dpV{p,e)\^^^ + EoCOs{e- 



-C0et-4>e)=^- 



(1) 



Our numerical analysis is based on the evolution equations, which describe the membrane potential y across the 
cellular membrane, along with a number of gating variables, collectively denoted as y, characterizing the con- 
ductance of various ionic channels. Symbolically, the system can be expressed as 
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dtV=-hon{V,y)/C+DV^V 
dtj = F{y,y), 



(2) 



where functions lion and Fare determined by different ionic currents 
in different models, C is the membrane capacitance, and D is the 
diffusion current coefficient. To demonstrate the results found in this 
paper are robust and essentially independent of precise ionic cur- 
rents, we use both Luo-Rudy modeP^ and Barkley modeP^. 

Based on equation (1), to describe the distribution of the mem- 
brane potential induced by CPEF near a circular obstacle, we apply 
CPEF at a weak strength in a two-dimensional quiescent medium. As 
shown in Fig. la (Luo-Rudy model) and Ic (Barkley model), the 
mechanism of WEH in response to CPEF rests on the induced de- 
polarization (red region) and hyper-polarization (blue region) near 
the obstacle. In both Luo-Rudy model and Barkley model, comparing 
to the dipole-like patterns induced by UEF^^'^^"^^ (Fig. lb and Id), the 
patterns induced by CPEF have two novel characters: one is that, they 
rotate synchronously with the rotating CPEF; the other is that, their 
patterns are similar as Chinese "ancient Taijitu". Depending on 
attaching to the obstacle or not, we divide both de-polarization 
and hyper-polarization into "Head" and "Tail" (Fig. la). 

The stably- rotating and fantastically- shaped de-polarization and 
hyper-polarization elucidated above can be used to unpin anchored 
spirals if the CPEF's strength Eq exceeds a certain value and angular 
frequency co^ is tuned to a proper value. In cardiac tissues, anchored 
spirals may have clockwise or counter-clockwise rotating directions. 
Therefore, to get a more comprehensive understanding about unpin- 
ning by CPEF, we discuss these two types of rotating directions 
separately. 

Firstly, we numerically simulate a clockwise rotating anchored 
spiral at the angular frequency oj^ and discuss its unpinning mech- 
anism by CPEF. In Luo-Rudy model as shown in Fig. 2a, at t = 0, is 
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Figure 1 | Distribution of the membrane potential induced by CPEF and 
UEF. (a), CPEF in Luo-Rudy model, Eq = 0.05 V/cm, = 0.2 rad/ms. 
(b), UEF in Luo-Rudy model, Eq = 0.05 V/cm. (c), CPEF in Barkley model, 
Eq = 0.05, We = 4. (d), UEF in Barkley model, Eq = 0.05. In this and 
following numerical simulations of Luo-Rudy model, the obstacle size R = 
0.32 cm, and the excitability is set as same as in Ref. 38. In this and 
following numerical simulations of Barkley model, the obstacle size R= 3y 
and unless otherwise specified, the excitability is set as a = 0.8, b = 0.06. 
The obstacles are all applied additional no-flux boundary conditions as 
equation ( 1 ) shows. The red dotted arrows represent the directions of 
electric field. The red curved arrows mean CPEFs rotate counter-clockwise. 
The red and blue regions around obstacles demonstrate de-polarizations 
and hyper-polarizations, respectively. The two black double-headed 
arrows indicate how we divide the "Head" and "Tail". The patterns in (b) 
and (d) which we get from simulations are the same as in Fig. la of Ref. 28. 



the initial phase of CPEF relative to x axis, (/>5 is the initial phase of the 
anchored spiral front relative to x axis. In cardiac tissues, (/>e is always 
given at a certain value, while (/>5 would be arbitrary. But in order to 
keep wave patterns simple and get a convenient structure analysis by 
numerical simulations, we choose (/>e is arbitrary but restrict (/>5 to 
zero. Furthermore, we can define the initial phase difference between 
(/>g and (/>5 as Acj) = (j)^ — to simply demonstrate the initial config- 
uration of CPEF and the anchored spiral. Hence, the configuration at 
t = 0 in Fig. 2a can be defined as a given A(/>. Then at t = 20 ms, a new 
wave N has been nucleated and begins to collide and merge with the 
anchored spiral S. Later at t = 40 ms, the colliding parts have 
detached from the obstacle, and form a new free spiral S'. 
Although there is another new wave N' nucleated by CPEF, its 
evolvement does not affect the final result. So after applying CPEF 
for a period in which an anchored spiral rotates one round (2n/Ws), 
there will be only S' left. And after ceasing CPEF for another period of 
In/ CDs, at ^ = 105 ms, S' does not re-pin to the obstacle but still keeps 
rotating freely. This is viewed as a "successful unpinning"^°. 

Additionally, the unpinning procedure above is much clearer in 
Barkley model, as shown in Fig. 2b. At the start, the configuration of 
CPEF and the anchored spiral S is introduced by a given A(/>. Then, at 
^ = 1 .0 when the phase of CPEF is 1 .OcOg + (/>e, a new wave N has been 
nucleated by the de-polarization. The one end of N, corresponding to 
the "Head" of de-polarization, is going to collide with S. Later att = 
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Figure 2 | Unpinning the clockwise rotating anchored spiral by CPEF. 

(a), In Luo-Rudy model, the angular frequency of spiral co^ = 0.136 rad/ 
ms. The CPEF's strength Eq = 0.7 V/cm and angular frequency cOg = 
0.1 rad/ms. CPEF is applied from t = Oto t = 46.2 ms. 0g is the initial 
phase of CPEF relative to x axis, and is the initial phase of the anchored 
spiral front relative to x axis and set as zero, (b), In Barkley model, the 
angular frequency of spiral co^ = 1.024. The CPEF's strength Eq = 1.8 and 
angular frequency co^ = 3.686. CPEF is applied from t = 0 to t = 6. N and 
N' represent different new waves nucleated by CPEF in different time. S 
and S' represent the initial anchord spiral and the new free spiral, 
respectively. White arrows are the propagation directions of waves. 
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(a) t=o 



t=8 ms 




Figure 3 | Unpinning the counter-clockwise rotating anchored spiral 
by CPEF. (a), In Luo-Rudy model, the angular frequency of spiral 
(jQs = 0.136 rad/ms. The CPEF's strength Eq = 0.7 V/cm and angular 
frequency cOg = 0.1 rad/ms. CPEF is appUed from t = 0 to t = 46.2 ms. 
(b), In Barkley model, the angular frequency of spiral co^ = 1.024. The 
CPEF's strength Eq = 1.8 and angular frequency cOg = 3.686. CPEF is 
applied from t = 0 to t = 6. 



1.8, the colliding parts merge with each other. And the other end of 
N, corresponding to the "Tail" of de-polarization, has detached from 
the obstacle, because its propagation along the boundary of obstacle 
is inhibited by both the "Head" of hyper-polarization and the refract- 
ory tail of S. Thereby, N and S form a new unpinned spiral S'. Because 
of the continuing CPEF, another new wave N' has been nucleated by 
CPEF, but N' has no effect to the final result. Finally, S' can be viewed 
as a successfully unpinned spiral. Therefore, we can recognize such 
can lead to successful unpinning, and call it the proper A(/>. 
In the next, we discuss the other mechanism about unpinning a 
counter-clockwise rotating anchored spiral by CPEF with a proper 
A(/>. In Luo-Rudy model (Fig. 3a), at the start of applying CPEF (t = 
0), an anchored spiral S is counter-clockwise rotating around the 
obstacle. Then at ^ = 8 ms, a new wave N is nucleated by CPEF. 
Later at t = 25 ms, S is unpinned from the obstacle. Finally, at t = 
105 ms, S rotates freely, which satisfy the requisite of successful 
unpinning. 

The clearer process can be seen in Barkley model (Fig. 3b). 
Similarly at the start, the configuration of CPEF and the anchored 
spiral is introduced by a proper Acj). Then at ^ = 2.0, a new wave N is 
nucleated due to the de-polarization induced by CPEF. Later ait = 
2.4, the anchored spiral S gradually falls into the "Head" of hyper- 
polarization induced by CPEF, which is at the opposite obstacle 
boundary of N. Because of the inhibition caused by the "Head" of 
hyper-polarization, S is unpinned. Furthermore, because of the 
inhibition caused by the "Tail" of hyper-polarization, S is driven 
further away from the obstacle. Although N is still rotating along 
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Figure 4 | The relations between /?u^ and cojcog in Barkley model. 

Different lines represent p^^ plotted against coJWs under different Eq. We 
choose interval [3.0, 4.0] for instance to reflect the existence of the optimal 
cojcos- Since the range centered around coJcOs = 3.6 has a relative large Puw 
even under a weak electric strength (e.g. Eq = 0.8, red dotted line), we 
adopt this ratio as the optimal coJcOs in following numerical simulations of 
Barkley model. 

the obstacle, it makes no effect to the final result and S will evolve to a 
successfully unpinned spiral. 

To summarize, we get two types of unpinning mechanisms by 
CPEF: for a given CPEF and excitability, with the proper A(/>, the 
rotating de-polarization and hyper-polarization induced by CPEF 
can lead to successful unpinning, corresponding to Figs. 2 and 3 
respectively. 

Therefore, we can consider A(/> is an important factor for successful 
unpinning. We define the whole range of A(/> which can lead to 
successful unpinning as the unpinning window {A(/>}unpin- Since A(/> 
is in the interval of [0, In), {A(/>}unpin can be normalized by 2n as 

p„„=MpxlOO%. (3) 

Pu^ also means the success rate of an arbitrary A(/> whether or not can 
lead to successful unpinning under the given CPEF and excitability. 

With given excitability which determines (Ds of an anchored spiral, 
Puw is highly related to the angular frequency co^ and strength Eq of 
CPEF. Among a certain range of co^, p^w can be optimized to its 
maximum. So we define the co^ giving maximal p^w as the optimal 
(De. Although different excitabilities (corresponding to different ojs) 
have different optimal co^ to optimize p^w^ the ratio of optimal co^ 
over given co^ keeps basically invariant. Thus we define this ratio as 
the optimal coJWs. We illustrate this relation between p^w and the 
optimal We/cOs, in the case of unpinning the counter-clockwise rotat- 
ing anchored spiral in Barkley model under a certain excitability and 
electric strength, as red dotted line in Fig. 4. Beside optimal ojJoJs, a 
proper electric strength Eq can also optimize Puw As also shown in 
Fig. 4, larger Eq makes larger p^^. And there is a limit beyond which 
increasing Eq makes no contribution to enhance Pu^v any more. On 
the other hand, a high electric strength harms hearts. So this limit 
value (e. g. Eq= 1.8 in Barkley model, as illustrated by black solid line 
in Fig. 4) could be used as the optimal electric strength. 

Now with the optimal settings of CPEF {Eq, cOg), we can examine 
the ability of CPEF to unpin anchored spirals in various excitabilities 
and compare it to the results by UEF. Firstly, to know how high the 
success rate (value of p^^v) is in various excitabilities, we take a part of 
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Figure 5 | A comparison of the success rate (value of /?u^) of CPEF and 
UEF in various excitabilities in Barkley model. The excitability is varied by 
changing b from 0 to 0.14 but fixing a to be 0.8. The blue soKd and red 
dashed lines represent the success rate of unpinning counter-clockwise and 
clockwise spirals by CPEF respectively, at Eq = 1.8, cOe = 3.6cOs. The 
black dotted line describes the success rate by UEF with the optimal 
strength Eq = 7, which we get from simulations and is the same as in Fig. 5a 
ofRef. 31. 

excitability region (a = 0.8, b = 0—0.14) for instance, as Fig. 5 shows. 
In this excitability region, p^w of CPEF reaches to an average of 80%, 
which is much higher than p^^ of UEF^\ It is clear that CPEF is very 
effective in successfully unpinning both counter-clockwise (blue 
solid line) and clockwise (red dashed line) rotating anchored spirals. 
Especially in weak excitabilities {b > 0.11), p^w of UEF is less than 
40%, but Puw of CPEF reaches 100%. This means that, regardless of 
A(^, CPEF can always successfully unpin anchored spirals in weak 
excitabilities. Moreover, in high excitabilities (^<0.06), UEF is failed 
at successfully unpinning any anchored spirals, but CPEF is still 
capable of successfully unpinning at an appreciable success rate 
(Puw > 40%). In addition, the optimal strength of CPEF (Eq = 1.8) 
required for successful unpinning is much smaller than the optimal 
strength of UEF (Eq = 7 in Ref. 31). 

Beside of high success rate, in the parameter space about excitabil- 
ity, the apphcation scope of CPEF, i.e. the extent where Puw > 0^ takes 
effect through the whole region where the medium exhibits excitable 
dynamics and supports spirals, and is much larger than the applica- 
tion scope of UEF^\ This is illustrated in Fig. 6: The area representing 
the application scope of CPEF (gray region) fulfills the whole region 
where spirals sustain (SW region), but the area representing the 
application scope of UEF (shaded region) is just within a part of 
SW region. Since cardiac tissues may distribute across SW region, 
CPEF is more applicable than UEF to unpin anchored spirals. 

In summary, we study the unique mechanism of WEH induced by 
CPEF, and find its outstanding ability to successfully unpin anchored 
spirals is better than UEF, at the higher success rate and larger 
application scope, even with a lower voltage. This is due to the unique 
characteristics of WEH induced by CPEF: in contrast to the dipole- 
like pattern induced by UEF^^'^^"^^, the pattern induced by CPEF is 
ancient-Taijitu-like (i.e. each of de-polarization and hyper-polariza- 
tion has "Head" and "Tail") and rotates synchronously with the 
rotating CPEF. Therefore, in the case of unpinning by de-polariza- 
tion, UEF would fail if the de-polarization is within the refractoriness 
of the anchored spiral and therefore cannot nucleate a new wave to 
unpin it. However, the rotating de-polarization induced by CPEF 
would escape from the refractoriness over time and nucleate a new 
wave easily. And the "Tail" of de-polarization induced by CPEF can 
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Figure 6 | A comparison of the application scope of CPEF and UEF in 
Barkley model. Three solid lines separate the whole parameter space about 
excitability into four regions marked by red double-headed arrows. In SW 
region, the medium exhibits excitable dynamics and supports spirals. For 
others, NW, RW and BI regions represent no wave, retracting waves and 
bi- stability, respectively. The shaded region represents successful 
unpinning by UEF with the optimal strength Eq = 7, which we get from 
simulations and is the same as in Fig. 2b of Ref. 31. The gray region 
including the shaded region represents successful unpinning by CPEF {Eq 
= 1.8, cOg = 3.6CO5). The vertical dotted line is corresponding to the chosen 
excitability used in Fig. 5. 

form an unpinned new spiral's tip naturally, since it is originally 
unpinned. In the case of unpinning by hyper-polarization, UEF 
would fail if the hyper-polarization is not within the attached region 
of the anchored spiral and therefore cannot unpin it. However, the 
rotating hyper-polarization induced by CPEF would meet the 
attached region of the anchored spiral over time and unpin it easily. 
Moreover, the "Tail" of hyper-polarization induced by CPEF can 
drive the unpinned spiral further away from the obstacle, to prevent 
re-pinning which is an important cause to unsuccessful unpinning by 
UEF. 

We believe, the outstanding ability of WEH induced by CPEF to 
unpin anchored spirals can be easily verified in cardiac tissues: 
instead of applying two DCs onto two pairs of field electrodes per- 
pendicular to each other to deliver UEF in the experimental prepara- 
tion of Fig. 5D in Ref. 16, one can apply two ACs onto these two pairs 
of field electrodes to realize CPEF in cardiac tissues, which is similar 
with the case in the Belousov-Zhabotinsky reaction of Ref. 35. We 
hope this lower-voltage higher- effectiveness approach may provide a 
better alternative to traditional therapies in terminating arrhythmia. 

Methods 

The Luo-Rudy modeP^ can be expressed as 

dt Cm (4) 

= Ino + Isi -\-Ik-\- Iki + hp + h, 

where y is the membrane potential, and lion is the total ionic currents which consist of 
a fast sodium current I^ay a slow inward current I^i, a time-dependent potassium 
current Ik, a time-independent potassium current Iki, a plateau potassium current 
iKp, and a time-independent background current 4. = 1 |iF/cm^ is the membrane 
capacitance, and D = 0.001 cm^/ms is the diffusion current coefficient. 

In the polar coordinate, equation (4) is integrated on a Np X Mq = 75 X 336 
uniform grid with no-flux boundary conditions via Euler method, and a five-point 
finite difference scheme is applied to compute the Laplacian term V^V. The space and 
time step are Ap = 0.02 cm, AO = n/16S and At = 0.001 ms. 
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The Barkley modeP^ can be expressed as 

du 



&~^u{l-u)[ u-^—\ +V^w 



(5) 



where u is the fast variable corresponding to the membrane potential, and v is a slow 
variable corresponding to the recovery process. The parameter £ determines the 
timescale of u which is fixed to 0.02. And the parameters a and h control the excit- 
ability of the medium: Larger a increases the action potential duration and larger h/a 
increases the excitation threshold. 

In the polar coordinate, equation (5) is integrated on a Np XMq = 91 X 396 
uniform grid with no-flux boundary conditions via Euler method, and a five-point 
finite difference scheme is applied to compute the Laplacian term Wu. The space and 
time step are Ap = 1/6, A6 = 7i/198 and At = 2 X 10"^ 
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